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ABSTRACT 

The chemical composition of high redshift galaxies is an important property which gives clues to their past 
history and future evolution and yet is difficult to measure with current techniques. In this paper we investigate 
new metallicity indicators, based upon the strengths of stellar photospheric features at rest-frame ultraviolet 
wavelengths. By combining the evolutionary spectral synthesis code Starburst99 with the output from the non- 
LTE model atmosphere code WM-basic, we have developed a code that can model the integrated ultraviolet 
stellar spectra of star-forming regions at metallicities between 1/20 and twice solar. We use our models to 
explore a number of spectral regions that are sensitive to metallicity and clean of other spectral features. The 
most promising metallicity indicator is an absorption feature between 1935 A and 2020 A, which arises from the 
blending of numerous Fe III transitions. We compare our model spectra to observations of two well studied high 
redshift star-forming galaxies, MS1512-cB58 (a Lyman break galaxy at z em = 2.7276), and Q1307-BM1163 
(a UV-bright galaxy at z em = 1.411). The profiles of the photospheric absorption features observed in these 
galaxies are well reproduced by the models. In addition, the metallicities inferred from their equivalent widths 
are in good agreement with previous determinations based on interstellar absorption and nebular emission lines. 
Our new technique appears to be a promising alternative, or complement, to established methods which have 
only a limited applicability at high redshifts. 

Subject headings: cosmology: observations — galaxies: high redshift — galaxies: evolution — galaxies: 
abundances — galaxies: starburst — stars: early-type 



1. INTRODUCTION 

It is now possible, using a variety of techniques, to assemble 
well-defined and large samples of galaxies at redshifts span- 
ning most of the Hubble time, from the present to z = 5. While 
many of the global properties of these galaxies — such as their 
luminosity function, clustering, and contribution to the star 
formation rate density — have been, or are, in the process of 
being characterised, there have been relatively few studies to 
date of their individual properties, such as their stellar pop- 
ulations and chemical enrichment, particularly at z > 1 . The 
simple reason is that, at these redshifts, the faintness of all but 
the brightest galaxies makes spectroscopic observations chal- 
lenging, even when using the largest aperture telescopes in 
the world in conjunction with highly efficient spectrographs. 
Nevertheless, a small number of detailed observations, which 
have targeted the brightest members of the population, have 
begun to provide some insight into the chemical properties of 

1 Based, in part, on data obtained at the W.M. Keck Observatory, which 
is operated as a scientific partnership among the California Institute of Tech- 
nology, the University of California, and NASA, and was made possible by 
the generous financial support of the W.M. Keck Foundation. 

- Current address: Isaac Newton Group, Apartado de Correos 321, 38700 
Santa Cruz de La Palma, Spain. 



high redshift galaxies. 

A number of authors, including Teplitz et al. (2000), Kobul- 
nicky & Koo (2000) and Pettini et al. (2001), have used neb- 
ular emission lines of [O II], [O III], and H/3 (the R23 method 
of Pagel et al. 1979) to measure the oxygen abundance of the 
ionized gas in galaxies at z — 3 selected with the Lyman break 
technique (Steidel et al. 1996). These studies have shown that 
the brightest of these Lyman break galaxies (LBGs) appear to 
be neither supersolar nor as metal poor as the damped Lya 
systems at comparable redshifts. Unfortunately, a degeneracy 
in the R23 calibration can lead to uncertainties of up to one 
order of magnitude in the derived oxygen abundance. More 
recently Shapley et al. (2004) have extended this work to us- 
ing [N II] and Ha (the N2 index discussed by Pettini & Pagel 
2004), to show that galaxies at z — 2 that are bright at rest- 
frame optical wavelengths are also metal-rich, with near-solar 
values of (O/H). However, one drawback of methods based 
on nebular emission lines from H 11 regions is that, at these 
redshifts, most of the lines used as abundance diagnostics fall 
at near-infrared wavelengths. In this spectral region obser- 
vations from the ground are made difficult by a multitude of 
strong night sky emission lines and by the limitations of cur- 
rent detector technology. 
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It is simpler, observationally, to perform optical rather than 
infrared spectroscopy and thereby study the rest-frame ultra- 
violet (UV) spectra of high redshift galaxies. Such spectra 
consist of the integrated light from the hot and luminous O 
and B stars in the galaxy (the stellar spectrum) on which are 
superimposed the resonant absorption lines produced by the 
interstellar gas (the interstellar spectrum). By exploiting the 
large wavelength coverage and high efficiency of the Echel- 
lette Spectrograph and Imager (ESI) on the Keck telescope, 
Pettini et al. (2002) obtained a high resolution and high signal- 
to-noise ratio (S/N) optical spectrum of the gravitationally 
lensed LBG MS1512-cB58 at z = 2.7276. From a detailed 
analysis of the interstellar absorption lines, these authors were 
able to determine, for the first time, the abundance pattern of 
many elements in a Lyman break galaxy. As well as indicating 
a relatively high metallicity, Z ~ 2/5Z Q , their results suggest 
a rapid timescale for the metal enrichment, of the order of a 
few hundred million years. Such in-depth studies are, how- 
ever, observationally very demanding and with present means 
can only be carried out on a few exceptionally bright galaxies. 

To overcome these difficulties we consider in this paper a 
novel technique for abundance determinations based on the 
rest-frame UV stellar spectrum. Such a technique has of 
course been applied extensively to the analysis of the optical 
spectra of old stellar populations, such as in elliptical galax- 
ies, for the last twenty years through the well-established sys- 
tem of 'Lick indices' (Burstein et al. 1984; Faber et al. 1985). 
However, its possible extension to the UV spectra of star- 
forming galaxies has not been addressed up to now. And yet, 
even at relatively low resolution and S/N, the rest-frame UV 
spectra of high redshift galaxies show a rich variety of stellar 
features which are mostly photospheric blends of absorption 
lines, as well as lines produced in stellar winds. It is thus 
worthwhile exploring to what extent the stellar spectrum can 
yield new metallicity probes (see, for example, the very recent 
work by Keel et al. 2004). 

As the integrated stellar spectrum from a whole galaxy has 
contributions from many different types of stars, it is not 
straightforward to extract information about the underlying 
stellar population. Nevertheless, progress in this direction can 
be made by employing spectral synthesis codes such as Star- 
burst99 (Leitherer et al. 1999, and references therein). Such 
codes consider a particular star formation law and stellar mass 
function and use theoretical evolutionary tracks to evolve the 
stellar population with time. Employing empirical libraries of 
stellar spectra to sum the contributions from individual stars, 
the codes can then predict the evolution of the population's 
integrated UV spectrum. By comparing an observed spec- 
trum with those synthesized with a range of different model 
parameters, it is then possible to constrain the properties of 
the underlying young stellar population [see Leitherer (2003) 
for a review] . 

The original empirical libraries adopted by Starburst99 
were compiled from observations of Galactic stars obtained 
by the International Ultraviolet Explorer (IUE) satellite. In 
order to assess the sensitivity of the stellar spectrum to metal- 
licity, a lower metallicity library was later appended following 
a Hubble Space Telescope (HST) programme to observe hot 
stars in the metal-poor Large and Small Magellanic Clouds 
(Leitherer et al. 2001). The success of this work demonstrated 
the potential of using spectral synthesis modelling to deter- 
mine the metallicity of stellar populations, as discussed in 
more detail below (§3.3). However, Galactic and Magellanic 
Clouds stars only sample a narrow range of metallicities. A 



major drawback for extending the empirical libraries of stel- 
lar spectra to other metallicities is that only relatively nearby 
stars can be observed with HST at the required spectral reso- 
lution and S/N. Even enlarging the libraries to include other 
Local Group galaxies would be challenging and observation- 
ally expensive while providing little increase in metallicity pa- 
rameter space. With no prospects in the foreseeable future for 
extending the empirical libraries to higher and lower metallic- 
ities, it may appear as if we have reached a limit to this type 
of analysis. 

Here we propose that a solution to this problem is to re- 
place the empirical libraries in Starburst99 with theoretical 
ones. A similar approach has been recently developed for the 
Lick indices by Thomas et al. (2003) and Thomas & Maraston 
(2003). Besides increasing the range of metallicities that can 
be investigated, theoretical spectra offer the additional ben- 
efits of a wider wavelength coverage and the possibility of 
altering the relative abundances of different elements, as well 
as the overall degree of metal enrichment, thus catering for a 
variety of star formation histories. This flexibility may be par- 
ticularly important when dealing with galaxies in the distant 
past, when the dominant mode of star formation may have 
been quite different from that commonly encountered in the 
nearby universe. Fortunately, the increasing sophistication of 
modern hot star model atmosphere codes makes it possible to 
predict the emergent UV spectra from stars of different struc- 
tural parameters and stellar wind properties. 

We have therefore used the non-LTE line-blanketed model 
atmosphere code WM-basic (Pauldrach et al. 2001), which in- 
cludes the effects of stellar winds and spherical atmospheric 
extension, to create theoretical libraries of stellar spectra over 
a range of metallicities (§2.1) and integrated these libraries 
into Starburst99 (§2.2). The outputs of Starburst99 mod- 
els thus produced compare favorably with more conventional 
models based on empirical libraries of stellar spectra (§2.3). 
The coupling of state-of-the-art codes from the fields of hot 
stars and starburst galaxies has allowed us, for the first time, to 
investigate the sensitivity of the integrated spectrum of a star- 
forming region to metallicity from Z = 0.05Z© to 2Z Q (§3). 
We consider in particular three blends of photospheric lines 
from OB stars— centered at 1370 A, 1425 A, and 1978 A— 
and conclude that the last one of these is probably the most 
suitable for stellar abundance determinations. In §4 we com- 
pare our models with Keck observations of two bright star- 
forming galaxies and find that the metallicities deduced from 
the stellar line indices agree to within a factor of ~ 2 with the 
values measured in the interstellar media of these galaxies us- 
ing established techniques. We discuss the results and draw 
our conclusions in §5. 

2. SPECTRAL MODELLING OF STAR-FORMING GALAXIES 

In this section we describe how we first generated theo- 
retical spectral libraries with WM-basic and then incorpo- 
rated them into Starburst99, in order to create a code capa- 
ble of modelling the UV spectral features of star-forming re- 
gions over a range in metallicities. Note that although WM- 
basic permits a free choice of metallicities, the evolutionary 
tracks and model atmospheres employed by Starburst99 are 
restricted to metallicities of 2, 1, 0.4, 0.2 and 0.05 Z Q . For this 
reason, we created theoretical libraries at these five metallici- 
ties. 

2.1. Generation of the WM-basic Spectral Libraries 
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2.1.1. The Non-LTE Line-blanketed Stellar Wind Code 
WM-basic 

In order to predict the emergent stellar flux at different 
wavelengths, it is necessary to produce a realistic model of 
the atmospheric structure of a star. For hot, massive stars 
this requires a substantial effort because the physical con- 
ditions in their atmospheres are complex and very different 
from those considered in standard models. They are dom- 
inated by the influence of the radiation field, which has an 
energy density larger than, or of the same order as, the en- 
ergy density of atmospheric matter. This has two important 
consequences. First, severe departures from local thermody- 
namic equilibrium (LTE) are induced, because radiative tran- 
sitions between ionic energy levels become much more im- 
portant than those caused by inelastic collisions between ions 
and electrons. As a result, we encounter non-LTE effects in 
the entire atmosphere of massive hot stars. Second, hydro- 
dynamic outflows of atmospheric matter (stellar winds) are 
initiated by line absorption of photons transferring outwardly 
directed momentum to the atmospheric plasma. The velocity 
fields of these outflows are mostly subsonic in the stellar pho- 
tospheres, where the weak photospheric lines are formed, but 
become highly supersonic above the photosphere, where the 
stronger spectral lines are formed. Both non-LTE and stel- 
lar winds affect the emergent energy distributions, ionizing 
fluxes and line spectra significantly and their effects need to 
be included in the spectral synthesis of hot stars (for a detailed 
discussion see, for example, Kudritzki 1998). 

Developing a code that solves the full complex set of cou- 
pled equations describing the physics of hot star atmospheres 
is a challenging numerical problem. These equations include 
the radiation hydrodynamics of the stellar wind (with the ra- 
diative line force from millions of spectral lines), the sta- 
tistical rate equations (for determining the occupation num- 
bers of different ionic species), the spherical radiative transfer 
equations for the spectral lines and the continuum, the energy 
equation and additional processes such as X-rays generated 
by shocks in the stellar wind outflow. It is not surprising that 
the development of such codes has become feasible only in 
recent years with increased computing power, improvements 
in numerical techniques, and expanded databases of atomic 
data. Examples of these sophisticated new codes are CM- 
FGEN (Hillier & Miller 1998), FASTWIND (Santolaya-Rey 
et al. 1997; Repolust et al. 2004), and WM-basic (Pauldrach 
etal. 1994, 2001). 

For the present work we chose to use the WM-basic code 
which has the important advantage of being extremely fast 
thanks to its use of a very advanced escape probability algo- 
rithm for the line transfer in expanding media. While this 
is an approximate treatment of the line transfer problem, it 
works very well when modelling the UV spectra of hot stars 
(Pauldrach et al. 2001). For the population synthesis work 
described here the speed of computation is crucial, since one 
needs to synthesize and add together the many different stel- 
lar spectral types that contribute to the integrated spectrum of 
a young stellar population. 

Before explaining in more detail how we used the code for 
our purposes, we note here that WM-basic treats the wind 
as homogeneous, stationary and spherically symmetric. Al- 
though winds are known to be time-dependent and unsta- 
ble, these assumptions seem to be reasonable approximations 
when modelling time- and population-averaged spectral prop- 
erties (see, e.g., Puis et al. 1993; Kudritzki & Puis 2000). WM- 
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FIG. 1. — Grid of WM-basic models. Models were computed far all five 
metallicities at each point marked in the Hertzsprung-Russell diagram. The 
grid points (filled circles) were chosen to lie on the solar metallicity evolu- 
tionary tracks of Meynet et al. (1994). The labels indicate the zero age masses 
considered while the lines show the tracks along which each star evolves. On 
the main sequence, stars of lower metallicity than solar are hotter; we there- 
fore considered additional points at higher temperatures (crosses) to sample 
this part of the Hertzsprung-Russell diagram. The other stellar parameters 
used when running the WM-basic models are summarised in Table 1 . 



basic also allows for the inclusion of the effects of the ion- 
izing emission from randomly distributed shocks in the stel- 
lar wind outflow using the approach developed by Feldmeier 
et al. (1997). Taking into account emission from shocks may 
be important for modelling correctly strong wind lines from 
high ionization stages such as O VI, O V, N V, and C IV (Paul- 
drach et al. 1994, 2001; Taresch et al. 1997), but photospheric 
and low ionization wind lines are unaffected by this process. 
Since in this first step of our work we will focus on the spec- 
tral synthesis of photospheric line features, we decided not to 
include the effects of shock emission at this stage. 

2.1.2. Computation of the Grid of Models 

WM-basic is distributed by A.W.A. Pauldrach as freeware 
software 3 on the world wide web and can be run on an IBM- 
compatible personal computer (PC) with a Microsoft Win- 
dows operating system, at least 350 MHz CPU, 64MB RAM 
and 200 MB free disk space. For the purposes of this work, 
it was installed on a 2 GHz Pentium PC in Cambridge and on 
a 1 GHz Athlon PC in Hawaii, both with 512MB RAM. De- 
spite the sophistication of the code, the operation of WM-basic 
is straightforward: the user specifies the stellar parameters of 
the star to be modelled and the program outputs the corre- 
sponding synthetic spectrum. We computed theoretical UV 
spectra for 87 grid points in the Hertzsprung-Russell (H-R) 
diagram at each of the five metallicities selected, with an indi- 
vidual model taking between 35 and 80 minutes to run. As we 
required a high spectral resolution, for improved accuracy, we 
requested output resolutions of 0.6 A and 0.25 A in the wave- 

3 http://www.usm.uni-muenchen.de/people/adi/Programs/Download.html 
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length ranges 900-1 150 A and 1150-2124 A respectively. In 
the remainder of this section we discuss the selection of the 
grid points, which are summarised in Table 1 and Figure 1, 
and the other input parameters used by WM-basic. 

As the purpose of our grid of models was to create li- 
braries of theoretical spectra from which Starburst99 could 
then synthesise the integrated UV spectra of star-forming re- 
gions, we first assessed which stars make a significant contri- 
bution to the total UV light. This is not immediately obvious 
because while hotter and more luminous stars produce prodi- 
gious amounts of UV flux they are much rarer and shorter- 
lived than lower mass stars. Simulations using Starburst99 
with its empirical libraries show that ~ 90% of the UV light at 
1450 A and 1800 A is generated by stars with zero age main 
sequence mass greater than 5M Q . For this reason, we re- 
stricted our grid to stars evolving from ones with masses of at 
least this value. 

To run realistic spectral synthesis models one must sample 
this upper part of the H-R diagram with sufficient grid points. 
As Starburst99 employs the evolutionary tracks of Meynet 
et al. (1994) to evolve the stellar population with time, we 
used these tracks as a guide for which grid points to select. 
We chose points corresponding to both zero age and evolved 
stars with initial masses of 120, 100, 80, 60, 50, 40, 30, 25, 20, 
15, 10 and 5 M Q , as illustrated in Figure 1. The evolutionary 
tracks were used to specify the full set of stellar parameters 
at each grid point, namely luminosity L, effective temperature 
r e ff and photospheric radius R. We list these values in Table 1 . 

It is also necessary to specify the stellar wind properties. 
For simplicity we first discuss stars of solar metallicity. The 
terminal velocity of the wind of each star was deter- 
mined from the empirical scaling formula of Kudritzki & Puis 
(2000) which relates Voo to the effective temperature and pho- 
tospheric escape velocity as follows: 

v OQ = C( r eff ) v esc , where ( 1 ) 

C(r eff ) = 2.65, r eff >21,000K 

c(r eff ) = 1.4, io,oook < r eff < 21,oook 
c(r eff )=i.o, r eff < io,oook. 

Note that the photospheric escape velocity is defined as 

v esc = [2^(l-D] - 5 (2) 

where g is the photospheric gravity, R is the photospheric ra- 
dius and r is the ratio of radiative Thomson to gravitational 
acceleration. This last parameter takes into account the effect 
of Thomson scattering in weakening the effect of the gravita- 
tional potential. 

We also had to estimate the mass loss rate M from each 
star in the grid. One of the predictions of radiation-driven 
wind theory [see Kudritzki (1998) for details] is that a wind 
momentum-luminosity relation exists such that 

logD 

mom = logD o +xlog(L/L ), (3) 

for 

D mom =Mv 00 (R*/R Q f\ 

where the coefficients Do and x vary with spectral type and 
(for O stars) with luminosity class. By using the empir- 
ical calibrations by Kudritzki & Puis (2000) and the revi- 
sions by Markova et al. (2004), we determined the wind mo- 
menta and hence mass-loss rates at most grid points. Note, 



however, that the wind momentum-luminosity relationship of 
eq. (3), which seems to hold for supergiants and O-type gi- 
ants, apparently breaks down in O-type dwarfs with luminosi- 
ties log(L/L ) < 5.2 . For these stars we adopted the rela- 
tion logD mom = 12 + 3 x log(L/L ) to account for their much 
weaker winds. The same formula was applied to B-type giants 
and dwarfs with luminosities log(L/L Q ) > 4.0 . For lower lu- 
minosities a very low value logD mom = 12 was adopted to sim- 
ulate a very weak, unobservable, stellar wind. Values of M 
and Voo for solar metallicity are listed in the last two columns 
of Table 1. 

Given that stellar winds are believed to be driven by the 
transfer of photon momentum due to the absorption and scat- 
tering of radiation by metal lines, it is expected that the wind 
properties will be sensitive to the metallicity. To test this as- 
sertion observationally, the properties of Galactic hot stars and 
those in the metal-poor Large and Small Magellanic Clouds 
have been compared (Puis et al. 1996). While there is clear 
evidence for smaller terminal velocities in SMC stars (Wal- 
born et al. 1995), further data are required to quantify the de- 
pendence of the mass-loss rate on metallicity. Nevertheless, 
there is an obvious trend for smaller wind momenta in lower 
metallicity O stars, and investigations are now underway to 
see whether this is also the case for B and A supergiants. In 
our models we adopted the best theoretical estimates for the 
scaling of and M with metallicity, as given by Leitherer 
et al. (1992), Vink et al. (2001), and Kudritzki (2002, 2003): 

/ z x0.13 

Voo = Voo, I — I and (4) 

M = M Q (-j . (5) 

Note that the hydrodynamic structure of the WM-basic 
models is calculated by accounting for a radiative line force 
in the equation of motion, parameterized by a 'line force mul- 
tiplier' . In order to calculate models with pre-specified values 
of the mass-loss rate and terminal velocity, it is necessary to 
estimate values of this parameter. To this end, we used the 
analytical approach by Kudritzki et al. (1989) in an iterative 
algorithm to provide the WM-basic code with force multipli- 
ers that yield the appropriate values of Voo and M as the result 
of the solution to the equation of motion. 

When running the grid at the five metallicities (2, 1, 0.4, 
0.2 and 0.05 Z Q ), we had the choice of using either (i) WM- 
basic' s solar reference abundances scaled by a constant fac- 
tor or (ii) our own choice of individual element abundances. 
We first consider the case of modelling star-forming regions 
within our Galaxy of 'solar' metallicity. Two possible sets 
of element abundances are the solar scale and the Orion neb- 
ula scale; we compare them in Table 2. The differences be- 
tween the two, particularly in some of the element ratios [e.g. 
(N/O), (Mg/O) and (Si/O)] are still a source of concern; pre- 
sumably they reflect unresolved systematic errors in either, or 
both, sets of measurements. In the present work, we adopted 
the Orion nebula values for elements included in the compi- 
lation by Esteban et al. (1998); column (3) in Table 2 lists the 
corresponding correction factors we applied to the default so- 
lar scale provided by WM-basic. For elements not included in 
the compilation by Esteban et al. (1998) we used the default 
solar scale in WM-basic. 

Of course there is no reason a priori why the mix of differ- 
ent elements should be the same as in the Orion nebula (or the 
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TABLE 2 
Orion Abundance Scale 



X 


12+log(X/H)ono„ a 


Correction b 


12+log(X/H) c 


Ref d 


[X/H]orion 


c 


8.49 


-0.029 


8.44 


1 


+0.05 


N 


7.78 


-0.139 


7.93 


2 


-0.15 


o 


8.72 


-0.112 


8.74 


2 


-0.02 


Ne 


7.89 


-0.189 


8.00 


2 


-0.11 


Mg 


7.43 


-0.150 


7.58 


3 


-0.15 


Si 


7.36 


-0.184 


7.56 


3 


-0.20 


S 


7.17 


-0.152 


7.20 


3 


-0.03 


CI 


5.33 


-0.175 


5.28 


3 


+0.05 


At 


6.49 


+0.093 


6.40 


3 


+0.09 


Fe 


7.48 


-0.025 


7.50 


3 


-0.02 



REFERENCES. — (1) Allende Prieto et al. (2002) — note that we adopt the C/O 
ratio from this source and apply it to our solar oxygen abundance; (2) Holweger (2001); 
(3) Grevesse & Sauval (1998). 



a Orion nebula abundances from the compilation by Esteban et al. (1998). 
b Logarifhmic correction factors applied to the internal WM-basic solar scale, in order 
to obtain the Orion abundances in column (1). 
c Solar abundance scale. 

^References for the solar abundance scale (see below). 
c [X/H] no„=log(X/H)-log(X/H) Q . 



Sun). Element ratios can (and do) vary as a function of metal- 
licity, reflecting the past history of star formation of a galaxy. 
Departures from solar ratios are plausible in high redshift star- 
forming galaxies, given that these objects generally support 
much higher star formation rates than our Galaxy did at any 
time in its past. Under these circumstances, elements syn- 
thesized by long-lived stars, primarily iron and other Fe-peak 
elements, may well be underabundant relative to the products 
of massive stars (oxygen and other a-capture elements), as 
found by Pettini et al. (2002) in MS1512-cB58. While one of 
the advantages of our fully theoretical approach to synthesiz- 
ing the UV spectra of star-forming galaxies is that we could in 
principle choose any mix of elements, for this first attempt we 
decided to simply scale the Orion nebula abundance pattern 
to different metallicities, keeping the relative element ratios 
fixed at their values in the Z = 1Z Q model listed in Table 2. If 
the results of the work presented here are sufficiently promis- 
ing, it will be possible in future to generate models with dif- 
ferent abundance ratios to, for example, determine the level of 
a-element enhancement in the galaxies under study. 

2.2. Integration of the Spectral Libraries into Starburst99 

Ours is the latest in a series of developments and improve- 
ments to expand the capability of the spectral synthesis code 
Starburst99 since it was first released in September 1998. The 
original version (see Leitherer et al. 1999) did allow some pa- 
rameters to be calculated at five metallicities (2, 1, 0.4, 0.2 
and 0.05 Z Q ) thanks to the inclusion of evolutionary tracks 
and model atmospheres at these metallicities. However, it was 
only possible to synthesize the emergent ultraviolet spectrum 
of a star-forming region of solar metallicity because empiri- 
cal stellar libraries were available only for Milky Way O-type 
stars. The inclusion of spectra of Galactic B stars in October 
1999 (de Mello et al. 2000) allowed the useful spectral range 
of the synthesized spectra to be extended to 1850 A where 
these cooler stars make a significant contribution to the inte- 
grated light. The next major development to the code was the 
introduction in December 2000 of a library of O-type spec- 
tra assembled from HST STIS observations of stars in the 
Large and Small Magellanic Clouds (Leitherer et al. 2001), 



giving the user the option to synthesize the UV spectrum of 
a population with a 'hybrid' metallicity of about 1/4 solar, al- 
beit over the restricted wavelength range 1250-1600 A. Most 
recently, Smith et al. (2002) created a grid of 230 expand- 
ing non-LTE line-blanketed model atmospheres, at 2, 1, 0.4, 
0.2 and 0.05 Z Q , and incorporated them into Starburst99 for 
the purpose of calculating realistic ionizing fluxes from young 
stellar populations. 

Smith et al. (2002) utilized two sets of model atmospheres: 
the WM-basic code by Pauldrach et al. (2001) for O stars and 
the CMFGEN code by Hillier & Miller (1998) and Hillier & 
Miller (1999) for Wolf-Ray et stars. Thus, their O-star mod- 
els were calculated with the same code used in the present 
study. Since Smith et al. (2002) focused on the emergent ion- 
izing continuum below 9 12 A and not on the line spectrum 
at longer wavelengths, the spectral resolution of their mod- 
els is not sufficient for the absorption line studies that are the 
subject of this paper. However, the continuous spectra pre- 
dicted by Smith et al. (2002) and by us are identical — unless 
Wolf-Rayet stars were important contributors to the contin- 
uum from the integrated stellar population, which is usually 
not the case. 

For the present work, we built a non-public version of Star- 
burst99 in which the empirical UV library of IUE and HST 
spectra was replaced by the grid of theoretical spectra gener- 
ated with WM-basic as described in §2.1. Most of the imple- 
mentation followed the same steps as the standard version of 
Starburst99, but with two significant differences, as we now 
explain. 

First, the theoretical spectra are not characterized by spec- 
tral types but by their effective temperatures and gravities. 
Therefore, we can immediately link them to the positions in 
the H-R diagram predicted by the stellar evolutionary tracks. 
This is a major advantage over empirical libraries because it 
avoids the intermediate step of having to adopt a spectral type 
vs. effective temperature calibration. Recent theoretical and 
empirical work (Crowther et al. 2002; Bianchi & Garcia 2002; 
Martins et al. 2002) has shown that the effective temperature 
scale of O stars needs to be revised downward from previous 
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FIG. 2. — Comparison of two Starburst99 + WM-basic spectra (red) with 'conventional' Starburst99 output spectra generated with empirical libraries of 
Galactic and Magellanic Clouds OB stars (black). The models shown are for continuous star formation (sampled at an age of 100 Myr when the UV spectral 
features have fully stabilized) with a Salpeter slope of the IMF between 1 and IOOMq and metallicities as indicated. The integrated spectra synthesized with 
empirical stellar libraries include narrow interstellar absorption lines which are not present in the fully theoretical spectra. The libraries of Magellanic Cloud stars 
used by Starburst99 only extend to 1 600 A; this accounts for the reduced range of the 'SB99-Empirical' spectrum in the lower panel. 



calibrations. Such a revision affects spectral synthesis mod- 
els based on the empirical libraries but not those generated 
theoretically. 

Second, model atmospheres predict both the idealized line- 
free continuum, and the apparent continuum which includes 
line blanketing. Therefore we can link the theoretical spec- 
tra to the evolutionary models without having to normalize 
them first. Empirical UV spectra are always affected by an 
unknown amount of the interstellar reddening. As a result, 
they must be rectified prior to adding them to a library. Then, 
the rectified spectra are recalibrated to flux units by fitting a 
theoretical low-resolution continuum at each grid point in the 
H-R diagram. This process introduces additional uncertain- 
ties which are avoided when a purely theoretical method is 
employed. 

We performed extensive tests with the synthesized spec- 
tra. As discussed below (§2.3), the computed photo- 
spheric features agree well with those generated with em- 
pirical stellar libraries. In addition, the stellar wind lines of 
C IV AA1548, 1550 and Si IV AA1393, 1402 show reasonable 
agreement. However, we believe that a more detailed study 
of all wind lines is warranted; since these lines are not uti- 
lized in the present analysis, their discussion is deferred to a 
future paper. Once work on the wind lines is completed, we 
will release an updated version of the Starburst99 code which 
includes the new theoretical libraries. 

With our new, fully theoretical, code we can now model 
the spectra of star-forming regions as a function of time, for 
either an instantaneous burst or a continuous star formation 
law and with a choice of the initial mass function (IMF) slope 
and mass range, for metallicities of 2, 1, 0.4, 0.2 or 0.05 Z Q . 

2.3. Comparison with the Starburst99 Empirical Libraries 

Before using our new code to investigate the effects of 
metallicity on the integrated stellar spectrum, we performed 
a 'sanity check' by comparing the outputs of two of our Star- 



burst99 + WM-basic models with the corresponding model 
spectra generated in the conventional manner, that is using 
empirical libraries of stellar UV spectra. Since the latter 
are available for two metallicities, corresponding to Galactic 
and Magellanic Cloud stars respectively, we compared them 
with the fully theoretical spectra generated with Z = 1 Z Q and 
0.2 Zq. 4 The comparison is illustrated in Figure 2 for our 
'standard' case of a continuous star formation model with 
a Salpeter IMF between 1 and 100 Mq (see §3). The spec- 
tra generated with the WM-basic libraries were smoothed and 
normalised according to the procedures described in §3.2. 

It is important to realise, when comparing the two sets of 
spectra, that they differ in their interstellar absorption. The 
empirical stellar spectra, obtained from observations of real 
stars, include the most prominent UV interstellar absorption 
lines which can be strong, particularly along sight-lines to the 
Large Magellanic Cloud. The WM-basic spectra, on the other 
hand, do not suffer from this contamination and are purely 
stellar. When this difference is taken into account, Figure 2 
shows that there is a reasonably good agreement between the 
spectra produced by the empirical and theoretical libraries at 
both metallicities. Some of the discrepancies, such as those at 
wavelengths greater than 1700 A in the solar metallicity mod- 
els, appear to arise from the continuum normalization rather 
than from intrinsic differences. Of particular interest for the 
present work are two blends of photospheric absorption lines, 
at 1360-1 380 A (which also includes the weak stellar wind 
line O V A1371) and 1415-1435 A, which Leitherer et al. 
(2001) showed to be sensitive to metallicity. Inspection of 
Figure 2 reveals a fairly good correspondence in these wave- 

4 The empirical libraries for Magellanic Cloud stars include stars from 
both the Large and Small Magellanic Clouds, giving a 'hybrid' metallicity 
Zmc ~ O.25Z0. We compare them with our closest theoretical libraries, for 
Z = 0.2Zq . In both the empirical and theoretical case, the evolutionary tracks 
used by Starburst99 are for Z = 0.2Zq . 
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TABLE 3 
Continuum Points 
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length regions. 

3. SYNTHETIC UV SPECTRA AT METALLICITIES Z = 0.05 - 2Z 

3.1. Computation of the Starburst99 Models 

Having confirmed that our new version of Starburst99 pro- 
duces results that are consistent with those generated by the 
empirical libraries, we can now turn to its application to the 
study of high redshift star-forming galaxies. In order to do 
so, it is necessary first of all to specify the mode of star for- 
mation. Starburst99 considers two limiting cases, either an 
instantaneous (i.e. unresolved in time) burst in which all the 
stars form at time t = and the stellar population then evolves 
passively, or continuous star formation at a constant rate. This 
distinction is important for our purposes because in the burst 
scenario the integrated spectrum changes rapidly with time, 
on timescales of a few Myr, reflecting the short lifetimes of 
the O and B type stars. In this limit, the emergent spectrum 
is more sensitive to the age than the metallicity of the stellar 
population. On the other hand, in the constant star formation 
limit, a quasi-equilibrium is reached within a few tens of Myr, 
after which there is little spectral variation with time (e.g. Lei- 
therer et al. 1999, 2001). Clearly, such a situation is much 
more favourable to the use of stellar features as metallicity 
indicators. 

Both modes of star formation considered by Starburst99 are 
idealised approximations to the real progress of star formation 
in galaxies, which is likely to be a series of bursts. How- 
ever, when analysing the integrated UV stellar output from 
a whole galaxy, the continuous star formation mode seems 
to be the better description, at least for most Lyman break 
galaxies, as we now explain. If a single burst followed by 
passive evolution were a valid approximation, we would ex- 
pect to see large numbers of post-starburst galaxies in which 
the most massive stars have already died. While still bright 
in the UV continuum thanks to the longer-lived B stars, such 
galaxies would lack the nebular emission lines whose strength 
is determined primarily by the hottest stars present. In other 
words, while the optical emission lines would fade away after 
~ 10 Myr, the galaxies would remain UV bright, and still meet 
the colour criteria tuned to the photometric selection of star- 
forming galaxies at redshifts z > 1 (Adelberger et al. 2004), 
for several tens of Myr. Galaxies with these characteristics 
seem to be the exception in the surveys by Steidel and collab- 
orators (Pettini et al. 2001; Erb et al. 2003). 

In any case, causality arguments alone argue for star for- 



mation timescales greater than 30 -40 Myr. This seems to be 
the characteristic dynamical timescale Td yn = 2r\jij(y of Ly- 
man break galaxies at z — 3, which typically have half-light 
radii rj/ 2 ~ \-2h~ 1 kpc and velocity dispersions a ~ 75 ± 
25 km s" 1 (Pettini et al. 2001). Similar values of Td yn apply to 
photometrically selected star-forming galaxies at z — 1 .5-2.5 
[the 'BX' and 'BM' galaxies surveyed by Steidel et al. (2004) 
and whose kinematics have been studied by Erb et al. (2003, 
2004)]. On the basis of these considerations, we only consider 
continuous star formation models in the present work. 

The second assumption which underlies all stellar popula- 
tion modelling is the shape and the mass limits of the IMF, 
which is usually approximated by a power law of the form 

dn 

<j>(m)= — =Cm- a , (6) 
dm 

(or a combination of such power laws), where n is the number 
of stars, m is their mass, C is a normalization constant and a is 
the IMF slope. The slope a = 2.35 at the upper mass end orig- 
inally proposed by Salpeter (1955) has stood the test of time 
(Kennicutt 1998) and seems to apply to high redshift galaxies 
too (Pettini et al. 2000; Steidel et al. 2004). Accordingly, in 
all the modelling presented here we adopted a universal IMF 
with a Salpeter slope between 1 and 100M Q (these being the 
usual limits used in Starburst99). However, in §4.3.1 we also 
consider the consequences of adopting different parameters 
for the IMF. 

3.2. Continuum Fitting and Resolution 

With the assumption of a continuous star formation mode, 
with a Salpeter IMF between 1 and 100M Q and a constant 
star formation rate, we used our new version of Starburst99 to 
compute synthetic spectra, at time intervals of 5 Myr from 5 to 
150 Myr, for each of the five metallicities considered. Before 
analysing these spectra, however, it is necessary to perform a 
couple of processing steps, as we now explain. 

In order to measure any spectral features (whether wind, 
photospheric, or interstellar lines) it is necessary to first nor- 
malize the spectra. This process removes the shape of the 
underlying continuum by fitting an appropriate function to re- 
gions deemed to be free of emission/absorption features and 
then dividing the spectrum by this function. Second, a com- 
parison between observed and synthesized spectra is best car- 
ried out at comparable spectral resolution so that, in general, 
it is necessary to smooth one spectrum to match the resolution 
of the other. The important point to be aware of is that contin- 
uum fitting and spectral resolution are intimately related. This 
is because in reality there are very few, if any, wavelength in- 
tervals in the UV spectrum of a young stellar population that 
are truly free of photospheric absorption lines. Thus, smooth- 
ing the spectrum invariably results in the depression of even 
relatively line-free regions and produces an apparent 'pseudo- 
continuum' which is lower than the real continuum level, now 
no longer discernible. 

The primary motivation of this work is to generate a set of 
fully theoretical spectra which can be compared with those of 
star-forming galaxies at high redshifts. The largest samples 
of such galaxies are from the spectroscopic surveys by Stei- 
del and collaborators (Steidel et al. 2003, 2004) conducted 
with the Low Resolution Imaging Spectrograph (LRIS) on 
the Keck telescopes. The resolution of this data set, FWHM ~ 
2.5 A in the rest frame, or ~ 500 km s" 1 at 1500 A, is the major 
source of line broadening, and dominates over the broadening 
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FIG. 3. — Examples of the continuum fitting process, showing the spline curves fitted to the 'pseudo-continuum' points listed in Table 3. The model spectra 
shown are for our standard case of continuous star formation (at age lOOMyr) with a Salpeter slope of the IMF between 1 and IOOMq and metallicities as 
indicated. 



from internal stellar processes, stellar rotation, and the veloc- 
ity dispersion of the stars within a galaxy. We therefore sim- 
ply convolved the synthetic spectra generated by Starburst99 
+ WM-basic with a Gaussian of FWHM = 2.5 A. In order to 
define the continuum level, we used both our model spectra 
and actual spectra of Lyman break galaxies to select 'pseudo- 
continuum' points which are listed in Table 3. For each model 
spectrum, we fitted a spline curve through the mean flux in 
each of these narrow windows and divided by this fit to pro- 
duce a normalized spectrum. Figure 3 shows some examples 
of the spline curves fitted to the models. 

3.3. The Sensitivity of the Stellar Spectrum to Metallicity 

Figure 4 presents the smoothed and normalised spectra of 
the 2, 1, 0.4, 0.2 and 0.05 Z models, lOOMyr after the on- 
set of star formation (by this time the stellar features in the 
composite UV spectrum have fully stabilised and show no 
residual time variation). Examination of the figure shows 
that the spectral appearance is strongly altered by the pres- 
ence of metals. The most obvious change, as the metallicity 
decreases, is in the declining strengths of both the absorption 
and emission components of the stellar wind lines, especially 
C IV AA1548, 1550. Such a marked dependence on metallic- 
ity has recently been reported in all the strongest stellar wind 
lines by Keel et al. (2004) in their study of star-forming re- 
gions in nearby galaxies at far-UV wavelengths with the Far- 
Ultraviolet Spectroscopic Explorer. While it may be tempt- 



ing to calibrate these features in terms of metallicity, unfor- 
tunately (for the present purposes) C IV AA1548, 1550 is nor- 
mally blended with strong interstellar absorption in the spec- 
tra of high redshift galaxies (see, for example, Fig. 4 of Pettini 
et al. 2003). Furthermore, the strengths and profiles of these 
lines respond not only to metallicity, but also to the slope and 
upper mass cutoff of the IMF (§4.3.1) and to time dependent 
dust obscuration (Leitherer et al. 2002). Although there have 
been attempts to interpret the C IV (and Si IV AA1394, 1403) 
lines in terms of metallicity (e.g. Mehlert et al. 2002), a proper 
treatment requires careful deblending of the different compo- 
nents (photospheric, wind, and interstellar) which contribute 
to these composite spectral features. These complicating fac- 
tors are much less important for the numerous blends of pho- 
tospheric absorption lines which make up most of the UV 
spectrum of star-forming regions and which can also be seen 
from Figure 4 to scale with metallicity. 

Following the inclusion of the low metallicity library of em- 
pirical stellar spectra into Starburst99, Leitherer et al. (2001) 
investigated the influence of metallicity upon the spectral ap- 
pearance of a star-forming region. They drew attention to two 
blends of lines near 1370 A and 1425 A (which they attributed 
to O V A1371 and Fe V AA1360- 1380, and to Si III A1417, 
C in A1427, and Fe V A 1430 respectively) and showed that the 
equivalent widths of these blends decrease by factors of 2-4 
as the metallicity of the stellar population drops from solar to 
~ 1 /4 solar. 
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FIG. 4. — Fully synthetic UV spectra of star-forming regions of metallicities 2-0.05 Zq. All the models shown are for our standard case of continuous star 
formation (at age 100 Myr) with a Salpeter slope of the IMF between 1 and 100 Mq. Each spectrum has been convolved with a Gaussian of FWHM= 2.5 A and 
divided by a spline fit to pseudo-continuum points (see Figure 3). 



We can verify and extend these results with our Starburst99 
+ WM-basic models. Figures 5 and 6 show the equivalent 
widths we measure for the '1370' and the '1425' line indices 
[integrated, as prescribed by Leitherer et al. (2001), over the 
wavelength ranges 1360-1380 A and 1415-1435 A respec- 
tively] for different metallicities and ages. These measure- 
ments confirm that the lines indices increase monotonically 
with metallicity over the range Z = 0.05 - 2Z Q considered 
here. Equally important is the fact that the indices are sta- 
ble after ~ 50 Myr from the onset of star formation, at which 
point a quasi-equilibrium is reached with approximately con- 
stant relative numbers of each type of star, even though indi- 
vidual stars continuously form and die. Dependence on metal- 
licity and stability with time are both necessary conditions to 
avoid a degeneracy between age and metallicity which would 
compromise the usefulness of a line index. 

Although our model spectra support the initial suggestion 
by Leitherer et al. (2001) that the 1370 and 1425 indices may 
be useful metallicity indicators, these two spectral features are 
far from ideal for this purpose. First, they are weak, with 
equivalent widths of only < 1.5 A at Z ~ Z Q ; therefore their 
use as metallicity diagnostics requires data of high S/N. A sec- 
ond complication is the fact that these lines are blends of tran- 
sitions from more than one element and therefore only give 
some (poorly defined) 'average' metallicity. It would not be 
possible, for example, to use them to investigate departures 
of different elements from their solar ratios. Finally, their 
far-UV wavelengths make them difficult to measure from the 
ground at redshifts lower than z ~ 1.7, thereby excluding the 
large numbers of star-forming galaxies now being discovered 



at these intermediate redshifts (Steidel et al. 2004; Abraham 
et al. 2004). 

3.4. The Region 193 5-2020 A 

Aware of these restrictions, we considered other UV stel- 
lar features which do not suffer from the same problems. The 
most suitable candidate we have found is in the range ^1900- 
2050 A. The broad blend of lines in this region was identified 
thirty years ago with the S2/68 Ultraviolet Sky-Survey Tele- 
scope on board the TD1 satellite which provided one of the 
first views of the sky at ultraviolet wavelengths. From spec- 
trophotometric measurements of early-type stars at a resolu- 
tion of 35 A, Thompson et al. (1974) reported the presence 
of a broad absorption region between 1900 A and 1940 A, 
which they attributed to line blocking by numerous photo- 
spheric Fe III transitions (Elst 1967). A more extensive study 
by Swings et al. (1976) confirmed that many of the absorp- 
tion features between 1800 A and 2200 A are indeed Fe III 
transitions. These authors also showed that the strength of 
the absorption is greater in B-stars of earlier types and higher 
luminosity classes. Nandy et al. (1981) observed a similar 
trend for B stars in the Large Magellanic Cloud and noted 
that the absorption is much weaker in O stars. In a conference 
proceedings, Morgan et al. (1986) addressed the question of 
whether the Fe III feature is sensitive to metallicity by study- 
ing stars in the metal-poor Small Magellanic Cloud. Their re- 
sults showed that, on average, the B star absorption is weaker 
at lower metallicities. 

In a study of the ultraviolet properties of nearby starbursts, 
Heckman et al. (1998) used IUE spectra of ~ 20 galaxies to 
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FIG. 5. — Variation of the 1370 line index with time and metallicity. 
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FIG. 7. — Variation of the 1978 line index with time and metallicity. 
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FIG. 6. — Variation of the 1425 line index with time and metallicity. 



create composite spectra with mean metallicities of ~ 0.16Z Q 
and ~ 1 . \Zq . Their work showed that not only is the Fe III ab- 
sorption clearly visible in the spectra of starburst galaxies but 
also confirmed that this feature is indeed stronger in higher 
metallicity environments. In addition to these encouraging 
preliminary indications, the absorption at ~ 1900 -2050 A 
overcomes many of the difficulties discussed above for the 
1370 and 1425 indices. It is broader and stronger, it should be 
mainly sensitive to the iron abundance, and is accessible from 
the ground for galaxies down to z ~ 1 . It therefore merits a 
closer look. 

Until now, it has not been possible to model this region of 
the spectrum, at any metallicity, because the empirical stellar 
libraries used by Starburst99 do not extend to these mid-UV 
wavelengths. This is not a problem, however, for WM-basic, 
and the output from our runs was deliberately chosen to reach 
beyond 2000 A (see Figure 4). The spectra shown in Fig- 
ure 4, which are the first to model the integrated spectrum of 
a star-forming region at mid-UV wavelengths, indeed confirm 
that the strength of the Fe III blend increases with increasing 
metallicity. 

To assess its variation with both time and metallicity we 
have defined a new index, based upon the equivalent width 
measured between two appropriate wavelengths. In choos- 
ing these wavelengths, we took into consideration not only 
our synthetic models but also the observed spectra of Lyman 
break galaxies, so as to avoid contamination with other spec- 



tral features. Interstellar absorption lines are fortunately not a 
problem, provided we avoid Al III AA1855, 1863 which can be 
strong (the much weaker Zn II+Mg I blend at A2026 is less of 
a concern). Potentially more serious is the nebular emission 
line C III] A 1909 which is clearly seen in the composite spec- 
trum of z — 3 galaxies generated by Shapley et al. (2003), al- 
though its strength can apparently vary substantially between 
different galaxies. 

For these reasons, we decided to restrict the definition of 
our new metallicity index to the equivalent width measured 
between 1935 A and 2020 A. This defines the '1978' index. 
In our synthetic spectra, these two integration limits appear 
to be at wavelengths where the stellar pseudo-continuum is 
recovered. Figure 7 shows how the index varies with time and 
metallicity. As was the case for the 1370 and 1425 indices, 
the 1978 index increases monotonically with metallicity and 
reaches a steady-state value after ~ 50 Myr, thus satisfying 
both requirements for a useful metallicity indicator. It can 
also be appreciated, by comparison of Figures 5-7, that at 
Z = Z Q EW(1978) is ~ 4 times greater than EW(1370) and 
EW(1425). 

4. COMPARISON WITH THE UV SPECTRA OF HIGH REDSHIFT 
GALAXIES 

To summarize our work so far, we have computed theo- 
retical spectra of OB stars using the stellar wind code WM- 
basic and integrated them into the spectral synthesis code 
Starburst99. With this combination, we have generated fully 
theoretical composite UV spectra of star-forming regions for 
a range of metallicities, from twice solar to 1/20 solar, under 
the simplest assumption of continuous star formation with a 
Salpeter IMF. We have analysed three blends of photospheric 
stellar lines, at 1370, 1425, and 1978 A, and found them 
to be potentially useful metallicity indicators, in that their 
equivalent widths increase monotonically with metallicity and 
change little with time after ~ 50 Myr (in the idealised ap- 
proximation of continuous star formation at a constant rate). 
The 1978 index, developed here for the first time, appears par- 
ticularly promising in that its equivalent width is larger than 
those of the other two, it is free of contaminating interstellar 
lines, and is thought to be due primarily to Fe lines. 

Of course, the validity of these photospheric UV indices 
must be tested against other established metallicity measures 
before they can be applied with confidence to the study of 
high redshift galaxies. At present such a comparison is only 
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TABLE 4 

Observations of the Rest-Frame UV Spectra of High Redshift Galaxies" 



Object 


Zsys 


n 


Instrument 


Grating/Grism 


Exposure time 


AA b s 


AA res t 


Resolution d 


SNR e 










(grooves mm - ') 


(s) 


(A) 


(A) 


(A) 




Q1307-BM1163 


1.411 f 


21.78 


LRIS B 


400 


5400 


3005- 7105 


1245 - 2945 


2.0 


37 


MS1512~cB58 


2.7276*5 


20.41 


ESI 


175 


16000 


3860-10500 


1040-2815 


0.37 h 


43 








LRIS R 


900 


11400 


4300- 6020 


1155-1615 


0.80 


75 








LRIS R 


1200 


3600 


5875 - 7185 


1580- 1925 


0.56 


24 



"All the observations were obtained with the Keck telescopes. 
b Systemic redshift. 

C B and R refer to the blue and red arms respectively of the LRIS spectrograph. 
d Spectral resolution at rest frame UV wavelengths. 
c Typical signal-to-noise ratio per resolution element. 
f From rest-frame optical nebular emission lines. 
g From the stellar photospheric C III A2297 line. 

h The resolving power of the ESI spectrum is R = A/AA = 5200, which corresponds to 0.37 A at a wavelength of 1900 A. 



possible in two cases, the galaxies MS1512-cB58 and Q 1307- 
BM1 163. These galaxies are unusually bright, either because 
they are intrinsically luminous (Q1307-BM1 163) or gravi- 
tationally lensed (MS1512-cB58), and have been observed 
over a range of wavelengths at high S/N, allowing measure- 
ments of element abundances mostly from interstellar absorp- 
tion and emission lines. While the preliminary indications 
from these tests are encouraging, as discussed below, clearly 
a larger sample of objects suitable for such a cross-check must 
be assembled for a comprehensive assessment of the method 
we have developed. We now consider each galaxy in turn; 
Table 4 gives details of observations of their rest-frame UV 
spectra with the Keck telescopes. The models shown in the 
comparisons are all for our standard case of continuous star 
formation with a Salpeter IMF at age 100 Myr. These models 
are available in full from Table 5. 

4.1. MS1512-cB58 

Thanks to gravitational lensing by the foreground cluster 
MS 15 12+36 at z = 0.373, MS1512-cB58 (or cB58 for short) 
is the brightest Lyman break galaxy known and thus one of the 
most extensively studied. Following its spectroscopic con- 
firmation as a distant galaxy by Yee et al. (1996) (the sys- 
temic redshift is z sys = 2.7276), Pettini et al. (2000) analyzed 
spectra obtained with the Low Resolution Imaging Spectrom- 
eter (LRIS) on the Keck telescope with a view to probing the 
galaxy's physical and chemical properties. A more detailed 
study of the kinematics and chemical abundances of the in- 
terstellar gas in cB58 was made possible by high resolution 
and high S/N ratio observations performed with ESI on Keck 
(Pettini et al. 2002). Here we compare both the LRIS and 
ESI data with the synthetic spectra generated by Starburst99 
+ WM-basic. 

Since the observations of cB58 are of unusually high res- 
olution (see Table 4), we first degraded them to the FWHM 
= 2.5 A of the synthetic spectra by convolution with Gaus- 
sian profiles of the appropriate widths. The spectra were 
then normalised by division by a spline fit to the continuum 
points listed in Table 3 [omitting two regions near 1275 A and 
1750 A where absorption lines from intervening absorption 
systems, unrelated to cB58, were identified by Pettini et al. 
(2000)]. We restrict our comparison between observed and 
synthesized spectra to the wavelength regions of the 1425 



and 1978 indices, and neglect the 1370 index which is com- 
promised by lines from two intervening absorption systems, 
Mg II AA2796,2804 at Zabs = 0.8287 and Si IV AA1394, 1403 
at z abs = 2.6606 (Pettini et al. 2000, 2002). 
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FIG. 8. — Comparison of the observed spectrum of the z = 2.7276 Lyman 
break galaxy MS1512-cB58 (black) with fully synthetic spectra produced by 
the Starburst99 + WM-basic combination (red) for five different metallicities, 
from twice solar to 1/20 of solar. The left-hand set of panels shows the region 
of the 1425 index (a blend of Si, C, and Fe photospheric lines — see text), 
while the right-hand set is for the Fe III blend near 1978 A. Each pair of 
panels is labelled with the metallicity of the synthetic spectrum shown. The 
y-axis is residual intensity (after normalization by the continuum level). 
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Figure 8 compares portions of the observed spectrum with 
the corresponding theoretical spectra computed for five metal- 
licities in the 1400-1450A (left) and 1898-2025A (right) 
wavelength intervals. We show the LRIS spectrum for the for- 
mer and the ESI one for the latter (although the ESI spectrum 
is in general the more useful of the two for such a comparison, 
it does suffer from poor flux calibration in the 1425 A region 
where two echellette orders overlap). The first conclusion to 
be drawn is that the Starburst99 + WM-basic combination 
does a remarkably good job at synthesizing ab initio the gen- 
eral character of these spectral regions, attesting to the high 
degree of sophistication achieved by these codes. Among the 
five metallicities considered, Z = 1Z Q and Z = 0.4Z Q are those 
which most closely match the observations in both the 1425 
and 1978 regions (giving comparable rms residuals between 
observed and model spectra). 

A metallicity between Z and 0.4Z Q is in good agree- 
ment with other measures. Specifically, the rest-frame opti- 
cal emission lines analyzed by Teplitz et al. (2000) with the 
Rn method of Pagel et al. (1979) imply an oxygen abundance 
in the ionized gas of 12 + log(0/H) = 8.39, or ~ 0.5 of our 
reference 'solar' value 12 + log(0/H) = 8.72 in Table 2 for 
the Orion nebula. Although less precise, arguments based on 
the strengths of the P-Cygni profiles of the C IV, Si IV, and 
N V wind lines also indicate a slightly sub-solar metallicity 
(Pettini et al. 2000; Leitherer et al. 2001). 

By far the most detailed assessment of element abundances 
in cB58 comes from the study of interstellar absorption lines 
from cool gas by Pettini et al. (2002). These authors found 
that among the elements covered in the ESI spectrum of 
cB58, those which are promptly released into the interstellar 
medium (Mg, Si, and S) all have abundances of ~ 0.4 solar 
with a scatter of about 0. 1 dex, comparable to the measure- 
ment error. On the other hand, the Fe-peak elements they 
observed have a mean metallicity Z = (O-II+qq^Zq, that is 
they are less abundant by a further factor of ~ 4, which Pet- 
tini et al. attribute to a combination of dust depletion and 
time delay in their release from Type la supernovae. Evi- 
dently, the iron abundance in the OB stars, as deduced here 
from the 1978 index, does not agree with the gas-phase abun- 
dance of interstellar Fe determined by Pettini et al. (2002), but 
is close to those of interstellar Mg, Si, O, and S. At face value, 
(Fe/H)oB stars ^ 4x (Fe/H) gas . This could be an indication that 
most of the underabundance of interstellar Fe in cB58 is due 
to dust depletion — an interpretation not favoured by Pettini 
et al., or that there is an offset in the calibration of the 1978 
index with metallicity, or a combination of both. An addi- 
tional contributing factor may be the unfortunate coincidence 
that, at the redshift of cB58, the 1935-2020 A wavelength re- 
gion is contaminated by telluric absorption; although this was 
divided out, we could only do so approximately and the pro- 
cess may have left some residual telluric absorption adding to 
the photospheric Fe III lines. 

In summary, we conclude that the analysis of photospheric 
UV lines from OB stars in cB58, using the methods devel- 
oped here, yields an estimate of metallicity which agrees with 
that of the interstellar medium to within a factor of ~ 2 (al- 
lowing for a moderate degree of dust depletion of the Fe-peak 
elements). 

4.2. Q1307-BM1163 

The combination of high luminosity and relatively low red- 
shift makes the TZ = 21.67, z= 1.411 galaxy Q1307-BM1 163 



TABLE 5 

Model Spectra for Continuous Star Formation with Salpeter IMF 
Wavelength Z = 0.05Z o Z = 0.2Z e Z = 0.4Z o Z = 1Z Z = 2Z 



1150.00 


1.025 


0.987 


0.951 


0.890 


0.825 


1150.18 


1.026 


0.988 


0.951 


0.889 


0.823 


1150.36 


1.027 


0.989 


0.951 


0.887 


0.821 


1150.54 


1.028 


0.989 


0.950 


0.886 


0.819 


1150.72 


1.028 


0.989 


0.949 


0.884 


0.817 


1150.91 


1.028 


0.988 


0.949 


0.884 


0.817 


1151.09 


1.028 


0.988 


0.949 


0.884 


0.818 


1151.27 


1.027 


0.988 


0.949 


0.885 


0.819 


1151.45 


1.027 


0.988 


0.949 


0.885 


0.820 


1151.63 


1.026 


0.987 


0.949 


0.886 


0.821 



NOTE. — The units of wavelength are A. The other five columns list values of 
residual intensity at five different metallicities, as indicated in the column headings. 
Table 5 is published in its entirety in the electronic version of the Astrophysical Jour- 
rial. A portion is shown here for guidance regarding its form and content. 
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FIG. 9. — Comparison of the LRIS-B spectrum of the z = 1.411 galaxy 
Q1307-BM1163 from the survey by Steidel et al. (2004) (black) with Star- 
burst99 + WM-basic synthetic spectra for five metallicities (red). Other de- 
tails are as in Figure 8. 



one of the brightest in the BM/BX sample of galaxies in the 
'redshift desert' (z — 1.4-2.5) recently assembled by Stei- 
del et al. (2004). These authors reported an oxygen abun- 
dance 12 + log(0/H) = 8.53 ±0.25 for the H II gas from the 
observed [N II] A6583/Ha emission line ratio (the N2 index 
of Pettini & Pagel 2004); this value corresponds to a metal- 
licity Z = (0.65^ 3)Z Q when compared with our reference 
12 + log(0/H) = 8.72. The profile of the C IV wind line is 
also consistent with near-solar metallicity (see Fig. 6 of Stei- 
del et al. 2004). 
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TABLE 6 

Summary of Abundance Measurements in MS1512-CB58 and Q1307-BM1 163 



Galaxy 


Redshift 


IS Abs. Lines'' 


Rn" 


N2 L 1425 index d 1978 index 6 


MS1512-CB58 


2.7276 


Z= (O.4±O.1)Z 
(Fem) =(0.11^) x(Fe/H) 8 


(O/H) = O.5(O/H) f 


Z = 0.7Z o (Fe/H)=O.7(Fe/H) 


Q1307-BM1163 


1.411 






(0/H)= (0.65^3) X (O/H) ■■■ (Fe/H) = 1.3(Fe/H) 



a Abundances from interstellar absorption lines (Pettini et al. 2002). 

b Oxygen abundance in the ionized gas (Teplitz et al. 2000), deduced from the R 21 index (R 2 3 = { [O II] + [O III] )/H/3) of Pagel et al. (1979). 
c Oxygen abundance in the ionized gas (Steidel et al. 2004), deduced from the N2 index (N2 = [N II] A6583/Ha) of Pettini & Pagel (2004). 
^Metallicity of OB stars, deduced from the equivalent width of the photospheric 1425 A blend [eq. (7)]. 
e Iron abundance of OB stars, deduced from the equivalent width of the photospheric 1978 A blend [eq. (8)]. 

^Note that the measurement of Teplitz et al. (2000) is scaled to the value of (O/H) adopted in the present work. Teplitz et al. (2000) do not provide an error estimate 
to the oxygen abundance. The typical accuracy of the R23 method is a factor of ~ 2 (e.g. Kennicutt et al. 2003). 
g This value is a lower limit, because an unknown fraction of the Fe-peak elements may be depleted onto interstellar grains. 



Figure 9 compares the LRIS-B spectrum of Steidel et al. 
(2004) with our five synthetic spectra, again in the 1400- 
1450 A (left) and 1898-2025 A (right) wavelength intervals. 
The 1360-1 380 A region is less useful in this case because, 
at the redshift of Q1307-BM1 163, it is close to the blue edge 
of the LRIS-B spectrum, where the S/N is low. As before, 
the observations were suitably smoothed to match the res- 
olution of the model spectra (in this case only a moderate 
amount of smoothing was required, from FWHM= 2 A to 
FWHM= 2.5 A) and normalised by reference to the contin- 
uum points defined in Table 3. Inspection of Figure 9 indi- 
cates that the synthetic spectrum with Z = 0.4Z Q is the clos- 
est match (lowest rms difference) to Q1307-BM1 163 in the 
1415-1435 A region, while the solar metallicity model repro- 
duces best the Fe III blend near 1978 A. The latter may be 
more reliable, given the higher S/N of the data at these wave- 
lengths. However, in either case the conclusion, as for cB58, 
is that the photospheric UV stellar indices give abundances 
that agree to within a factor of ~ 2 with those deduced from 
other indicators based on stellar wind lines and nebular emis- 
sion lines from H II regions. 

4.3. The Use of Line Indices 

When we assessed the outputs of our code against the spec- 
tra of cB58 and Q1307-BM1 163 in the preceding sections, 
we did so by identifying the synthetic spectrum, among five 
options in metallicity, that gives the smallest rms deviations 
from the observations. For a more straightforward determi- 
nation of metallicity, we give here functional forms to the re- 
lations between the equivalent widths of the 1425 and 1978 
line blends and log(Z/Z ). These were deduced by fitting the 
values measured from our synthetic spectra at the five metal- 
licities for which they were generated, as shown in Figure 10. 
The values of EW fitted are those from Figures 6 and 7, at 
time lOOMyr. We find: 



log(Z/Z ) = A x EW(1425) + B, where 



(7) 



A = 1.74 
A= 1.14 



and 
and 



B = -1.75 
B = -1.39 



EW(1425)<0.6A 
EW(1425)^0.6A 



and 



C = 0.33 and 
C = 0.42 and 



D = -1.94 
D = -2.50 



EW(1978K5.9A 
EW(1978) > 5.9A. 



As can be seen from Figure 10, both indices increase lin- 
early, or nearly so, with logZ over the factor of 40 in Z con- 
sidered in this work. It is important to realize that the re- 
lations given in (7) and (8) above are valid only for spectra 
of 2.5 A resolution normalised through the continuum points 
listed in Table 3. Given the broad, shallow nature of these 
photospheric blends, their equivalent widths are very sensi- 
tive to the continuum normalization which, as explained in 
§3.2, in turn depends on the spectral resolution of the data. 

Applying eqs. (7) and (8) to cB58, we find logZ = -0.16 
(Z = O.69Z ) from the measured EW(1425)=1.08 A, and 
logZ = -0.17 (Z = O.68Z ) from EW(1978)=5.37A. Although 
in good agreement with each other, these values are some- 
what higher than Z = (0.4±0.1)Z Q determined for the in- 
terstellar gas (§4.1). In Q1307-BM1 163 we can only mea- 
sure the 1978 index, because the 1425 blend is affected by 
a noise spike — see Figure 9. (This is one of the limita- 
tions in reducing a whole spectral feature to a single number. 
While admittedly more objective, this approach is however 
more vulnerable to noise than the spectral match discussed 
in §4.1 and §4.2. With the latter it is still possible to re- 
cover useful information from 'clean' portions of the spec- 
trum). We deduce logZ = +0.10 (Z = 1.3Z ) from the mea- 
sured EW(1978)=6.2A in Q1307-BM1 163, which again is 




z/z 



Z/Ze 



log(Z/Z ) = CxEW(1978) + D, where 



(8) 



FIG. 10.— EW-Z relation for the 1425 and 1978 indices. The dots are 
the values of equivalent width measured from the synthetic spectra at age 
lOOMyr, while the lines show the functional fits given by eqs. (7) and (8). 
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TABLE 7 
IMF Parameters 



Model 


Star Formation Mode 


IMF slope" 


M lip (M )» 


A 


Continuous 


2.35 


100 


B 


Continuous 


3.30 


100 


C 


Continuous 


2.35 


30 



a The parameter a in eq. (6). 
^Upper mass limit 



somewhat higher than Z = (0.65^5 3) Zq deduced from the 
[N II]/Ha emission line ratio from H II gas (§4.2). 

All the above abundance measurements in cB58 and 
Q1307-BM1 163 are summarised in Table 6. Considering that 
each of the established methods (interstellar absorption lines, 
R23, and N2) give values of metallicity which are only accu- 
rate to within a factor of ~ 2, the photospheric line indices 
developed here appear to be a promising alternative. On the 
basis of this very limited comparison, it would seem that the 
hot star models tend to slightly underestimate the strengths 
of the photospheric blends at 1425 A and 1978 A at a given 
metallicity. It would be interesting, in future, to verify if this 
is indeed the case with a larger set of galaxies where both 
stellar and H II region abundances have been measured. 

4.3.1. Sensitivity of the Line Indices to the IMF parameters 

Finally, we briefly consider the sensitivity of the 1425 and 
1978 indices to the IMF parameters. To this end, we gen- 
erated a new set of synthetic spectra as described in §3, but 
with different IMF prescriptions, and remeasured values of 
EW(1425) and EW(1978) at the five metallicities considered 
in this work. Table 7 gives details of the IMF parameters con- 
sidered. Model 'A' is our standard case of continuous star 
formation with a Salpeter IMF at age lOOMyr. Model 'B' has 
a steeper power-law slope {a = 3.30 in eq. 6), and is thus close 
to the Miller & Scalo (1979) IMF at the high mass end. Model 
'C has the standard Salpeter slope but the IMF is truncated at 
M up = 30M Q . Thus, both models B and C have comparatively 
fewer high mass stars than model A. 

As can be seen from Figure 1 1 , decreasing the proportion of 
massive stars consistently increases the equivalent widths of 
the two blends (at a given metallicity), reflecting the fact that 
these spectral features are strongest in late O and early B-type 
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FIG. 1 1 . — Dependence of the 1425 and 1978 indices on the IMF param- 
eters. Model A is for our standard case of continuous star formation with 
a Salpeter IMF at age 100 Myr. Models B and C have comparatively fewer 
massive stars than model A (see Table 7). The lines here simply join the dots, 
to guide the eye. 
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FIG. 12. — Normalised, synthetic UV spectra at solar metallicity for the 
three sets of IMF parameters listed in Table 7. Model A is for our standard 
case of continuous star formation with a Salpeter IMF at age 100 Myr. Mod- 
els B and C have comparatively fewer massive stars than model A. 



stars. However, there are other spectral changes which accom- 
pany these IMF alterations. The most obvious one (see Fig- 
ure 12) is in the character of the C IV AA1548, 1550 line which 
in models B and C no longer exhibits a discernible P-Cygni 
component, as already noted by Pettini et al. (2000) and Stei- 
del et al. (2004). The empirical observation that Lyman break 
galaxies normally show such P-Cygni features — as demon- 
strated by the composite spectrum of Shapley et al. (2003) — 
argues against significant departures from a Salpeter IMF (at 
least for masses M > 5M Q ) even at these early epochs. 

5. DISCUSSION AND CONCLUSIONS 

Rest-frame UV spectra of galaxies at redshifts z > 1 are now 
being routinely obtained with large ground-based telescopes. 
Current samples already include thousands of galaxies with 
spectroscopic redshifts and their numbers will continue to 
grow in the years ahead. Encoded in all these spectra is in- 
formation on the chemical composition of the galaxies which 
is one of the key physical properties for understanding their 
previous history and future evolution. Potentially, knowledge 
of the relative abundances of the chemical elements in these 
galaxies can tell us about the accumulation of the products 
of previous stellar generations, the consumption of gas into 
stars, the pace of star formation activity, the initial mass func- 
tion, and the importance of inflows and outflows in shaping 
the chemical evolution of galaxies in the distant past. 

In this paper we have taken the first steps towards decod- 
ing some of this information. In order to do so, we have 
found it necessary to synthesize ab initio the spectra of stars 
at the upper end of the H-R diagram and use them, via a spec- 
tral synthesis code, to predict the emergent UV spectrum of a 
composite stellar population as a function of metallicity. This 
approach was dictated by the practical difficulties with cur- 
rent facilities to assemble such a library of stellar spectra di- 
rectly from observations and was made possible by the high 
degree of sophistication achieved by codes such as WM-basic 
and Starburst99. A further advantage of a fully theoretical 
method such as ours is that one has complete control over the 
wavelength coverage and resolution of the synthesized spec- 
trum, as well as over the metallicity and relative abundances 
of individual elements. 
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The results of this initial study are encouraging, in that the- 
oretical and empirical spectra resemble each other closely and 
in detail — while the models may not be perfect yet, they pro- 
vide a remarkably good representation of reality. 

We have focussed our attention on three blends of photo- 
spheric stellar lines, near 1370, 1425, and 1978 A, and found 
that they are indeed mostly sensitive to metallicity. In the 
idealized case of a continuous star formation activity at a con- 
stant rate, the equivalent widths of these three features vary 
little with time after approximately 50Myr, and increase by 
factors of <~ 4-6 as the metallicity of the stars increases from 
1/20 of solar to twice solar. Of the three blends, the one at 
1978 A may be the most useful because: (a) it is the strongest, 
with an equivalent width EW(1978)= 5.9 A at solar metallic- 
ity; (b) the wavelength interval 1935-2020 A over which it is 
defined is free of strong interstellar lines; (c) its rest-frame 
wavelength makes it accessible to ground-based spectroscopy 
over the redshift range z = 1 — 3 where recent surveys for high 
redshift galaxies have been prolific; and (d) it is due primar- 
ily to Fe III lines and thus measures the abundance of a single 
element rather than a combination of several. 

We have only been able to test the validity of our method 
in two cases, by comparing the synthetic spectra we gener- 
ated with those of two bright high redshift galaxies, MS 15 12- 
cB58 at z = 2.7276 and Q1307-BM1163 at z = 1.411. We 
find that in these two cases the metallicities deduced from the 
photospheric lines agree, to within a factor of <~ 2, with those 
derived from analyses of interstellar absorption and emission 
lines. More extensive comparisons will allow us to refine the 
calibrations of our line indices in terms of metallicity. 

The principal conclusion of this work is that photospheric 
UV lines offer a much needed complementary avenue to de- 
termining the degree of metal enrichment in high redshift 
galaxies. The main weaknesses of the method we have de- 
veloped are: (a) It is based on spectral features that are broad 
and shallow. Thus, while they can be detected even in spec- 
tra of modest resolution, their measurement does require a 
moderately high signal-to-noise ratio, S/N> 15 per pixel at 
R ~ 1000. (b) The equivalent widths of the features in ques- 
tion depend on the mix of spectral types present in a star- 
forming region. Thus, their strengths would vary with time 
in isolated bursts of star formation and, to a lesser extent, 
with the initial mass function. However, to date, neither single 
bursts nor departures from a universal IMF (at the upper mass 
end) have been found to be common in high redshift galaxies. 

The main strength of the technique is its complementarity 
to other abundance measures: it is applicable precisely where 
other techniques fall short. Abundance determinations from 
interstellar absorption lines require a combination of spectral 
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resolution and S/N that is not achievable with current means 
for the vast majority of high redshift galaxies. Circumventing 
this problem by co-adding spectra of many galaxies is fraught 
with dangers (e.g. Savaglio et al. 2004), since the absorption 
lines are normally on the flat part of the curve of growth, 
where their equivalent widths are much more sensitive to the 
velocity dispersion of the gas than to its column density (and 
therefore metallicity). Furthermore, most abundant elements 
are depleted onto dust in the interstellar medium; correcting 
for unknown and variable amounts of dust depletion is a sig- 
nificant complication which has not yet been fully resolved 
even in QSO absorption line studies (e.g. Vladilo 2004). 

Emission lines from H II regions, on the other hand, when 
observed at high redshifts are limited to the strongest lines 
which sample only a few elements (primarily oxygen). The 
so-called 'strong line methods' give estimates of the oxygen 
abundance that are only accurate to within a factor of ~ 2 [see, 
for example, the discussion by Pettini & Pagel (2004)], com- 
parable to the accuracy of the photospheric UV line method 
explored here, at least on the basis of initial tests. Further- 
more, there are only a few redshifts intervals at z > 1 where 
the nebular lines required for abundance measurements all fall 
within near-IR windows accessible from the ground, and even 
then contamination by strong sky emission lines is often a lim- 
iting factor to the measurement of reliable line ratios. 

Looking ahead, there are thus strong motivations for con- 
tinuing the work we have begun by: (a) incorporating fu- 
ture improvements in WM-basic into our hybrid code and, 
(b) most importantly, performing more extensive cross-checks 
between abundances deduced from nebular emission lines and 
with our photospheric UV absorption line technique, in cases 
where both methods are applicable. With a reliable calibration 
of the latter in hand, we would then have the means to mea- 
sure stellar abundances in a wholesale manner over a wide 
range of cosmic epochs. 
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Spectral Modelling of Star-forming Regions in the UV 



TABLE 1 
Grid of WM-Basic Models 
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log (L/L ) 


7k 


log (g/cm s 2 ) 
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(M ) 




(K) 




(*©) 


(M yr 1 ) 


(kms- 1 ) 


119.98 


6.25 
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4.333 


12.39 


0. 1 8E-04 


4076 


119.61 
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0.17E-04 
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NOTE — The wind parameters M and Voo refer to the values used for models at solar metallicities. For the remaining metallicities, these parameters were scaled 
according to eqs. (4) and (5). 



